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ABSTRACT 

Context. Oscillations in rotating stars with frequency a- of the same order or smaller than the rotation rate CI cannot be described by a single 
spherical harmonic due to the effect of the Coriolis force. This is a serious complication which is usually treated by writing the eigenfunctions 
as a (truncated) sum of different spherical harmonic degrees for a given m- value, or by neglecting the 9 part of the rotation vector (Traditional 
Approximation) . 

Aims. We aim for a more adequate treatment of the coupling with higher angular degrees for low frequency oscillations of rotating stars by 
taking the Coriolis force fully into account, so that the coupling is included in the analysis (up to the grid resolution) and to compare the results 
with analyses based on the above mentioned approximations. 

Methods. To this end a new, more efficient version of a 2D(r,0) implicit oscillation code was developed in which no a priori assumptions about 
the 6 variation of the eigenfunctions is made, enabling a better treatment of the rotational truncation of g-modes near the stellar poles, and for 
which ~ 150 gridpoints in 9 are feasible. We test the code by comparing the simulated oscillation spectrum with that obtained by the MOST 
photometry for the ^ 6 Mq Be star HD 163868. 

Results. We find both prograde and retrograde overstable modes (although more prograde than retrograde modes) and confirm the existence of 
low degree odd r-modes, destabilised by the /c-mechanism. The ultra-low frequency modes that could not be explained in a previous analysis 
are interpretated as (high degree) retrograde m=l modes with \a\ ~ CI i.e. cr ~ 0. A reasonably good fit to the observed oscillation spectrum is 
possible if we assume that only even modes are observed (no unstable r-modes visible). This requires a nearly equator-on view of the observed 
star, consistent with the measured high v sin ; value of 250 km/s. 
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1. Introduction 



In a previous paper d Savoniie 2005 ) the pulsational stability 
of uniformly rotating main sequence stars in the mass range 
3-8 Mq was studied by applying the traditional approxi- 
mation (TA). Applying the TA renders the oscillation equa- 
tions fully s eparab le, greatly simplifying the solution, see e.g. 



Unno et al. (Il989l) . The earlier calculations indicated that the 



/(•-mechanism in rapidly rotating stars can not only destabilise 
gravity (g) modes but, more surprising perhaps, also rotational 
r-modes. The TA analysis shows that the lowest degree odd 
(regarding the symmetry of the density and pressure perturba- 
tion with respect to the stellar equator) r-modes with low m- 
values have relatively large eigenvalues A for the angular part of 
the eigenfunctions. This indicates relatively large buoyancy ef- 
fects for these r-modes, almost comparable to that of g-modes. 
The compressions and expansions in the driving zone (near 
r ^ 2 X 10^ K) appeared indeed sufficiently strong for destabi- 
lization of some of these modes by the y-m echanism. This re- 
sult was corroborated by the (TA) analysis of lTownsendl (I2005I) 
and by the truncated spherical harmonics expansion method by 
Lee ( 20061) . Unfortunately the simplified TA analysis cannot 



take the coupling between different /-components through the 
Coriolis force adequately into account. 



Savonije & Papaloizoul (119970 developed an implicit two 



dimensional (2D) oscillation code, which solved the linear 
forced oscillation equations including the full Coriolis force 
(but neglecting the distorting centrifugal force), and applied it 
to the tidal evolution of massive binary systems. They found 
resonances with g-modes and r-modes and indications of inter- 
actions between the oscillations in the core and envelope when 
the forcing frequency cr (in corotating frame) was in the inertial 
regime (|a"| < 2f2), where Q. is the stellar rotation rate. 

Here we apply a similar implicit 2D code to the pulsational 
stability analysis of an evolved {Xc = 0.4) 6 Mq main sequence 
star rotating uniformly with O/Qc = 0.6, where Q^. is the sur- 
face break-up speed. By applying a tidal force on the star with 
varying (complex) driving frequency one can search for res- 
onances with free oscillation modes and determine their pul- 
sational stability. By varying the (sign of) the imaginary part 
of the forcing frequency, and searching for a maximum re- 
sponse, we can determine the damping- or growth rate of the 
mode. This procedure is different from the usual stability anal- 
ysis whereby one numerically integrates a (one dimensional) 
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eigenvalue problem by iterating until one reaches a consistent 
solution. Explicit integration like this, whereby the boundary 
conditions are only fulfilled after a number of iterations, is 
likely to be impractical for the set of 2D partial differential 
equations describing stellar oscillations with the full Coriolis 
force. Therefore we apply a 2D code which numerically inte- 
grates the oscillation equations implicitly. The boundary con- 
ditions are applied along the rotation axis, along the stellar sur- 
face and back along the stellar equator Starting at the stellar 
surface on the rotations axis, using the local boundary con- 
ditions together with the difference version of the oscillation 
equations, we start the elimination procedure that continues 
along all the points of the (r, 0) grid towards the stellar cen- 
tre. By finally applying the central boundary conditions we can 
then solve the unkown perturbations there and start a backsub- 
sitution in reverse order towards the stellar surface. The thus 
calculated result yields for each adopted forcing frequency (in 
the resonance search) a solution consistent with all boundary 
conditions, but it represents in general not a free oscillation 
mode. However, when the frequency gets close to resonance, 
the forcing term becomes insignificant and the forced solution 
approaches that of the free oscillation mode. The implicit na- 
ture of this 2D integration causes the calculations to be memory 
and cpu-time expensive as the calculations involve many ma- 
trix inversions during the elimination procedure, the more so 
if one requires high grid resolution to allow adequate coupling 
with higher degree, short wavelength components. In order to 
improve the speed, and thus allow higher resolution of both the 
radial and angular part of the oscillatory response, a new ef- 
ficient implicit 2D-code was developed based on a three-level 
difference scheme, see below. 

We compare our numerical results with the MOST obser- 
vations of the Be star HD 163868 and the calculated oscilla- 



tion spectrum by Walker et al.l (12005.) and iDziembowski et al 
(l2007h . 

2. Non-adiabatic Oscillation Equations 

We introduce spherical coordinates (r, 6, (p) centered at the stel- 
lar centre, whereby = corresponds to the stellar rotation 
axis, and solve the non-adiabatic linear oscillation equations 
with (the real part) of the spherical harmonic forcing potential 



i((7t-mip) 



= A r' P'^(0) e 



i(iTt-m(^) 



(1) 



where A is a constant and P"(0) the associated Legendre func- 
tion. We define m to be always positive and cr the o scillation 
freque ncy in the inertial frame, see e.g. eqs. (l)-(4) in lSavoniie 
( l2005b . To keep the problem tractable we ignore the centrifu- 
gal distortion of the (rapidly) rotating star but include the full 
Coriolis force on the oscillations. Because the unperturbed 
state of the star is assumed spherically symmetric and we look 
for a stationary solution, the time and (f part of the perturbed 
solution is separated as oc q^^'^-'^'p) and taken out of the equa- 
tions below. For better efficiency we pre-eliminate the if com- 
ponent of the displacement vector ^ together with the perturbed 
radiative flux components F'g and F'^. Hence only the four 
unknown perturbations ^e, P' IP and p'/p remain, where 
the symbols have their usual meaning and a prime denotes the 



Eulerian perturbation of a quantity. We ignore perturbations of 
the convective energy flux. By defining & = cr - mQ., the os- 
cillation frequency in the corotating frame (whereby cr < 
corresponds to retrograde oscillations and o- > to prograde 
oscillations in stellar frame), the equation for the perturbed ra- 
dial motion of a mass element can be expressed as 
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while the perturbed equation for mass conservation follows as 
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The perturbed energy equation with radiative diffusion can 
be written as (with i the imaginary unit) 

(Fj - 1) 



P \ p / \ dr dr 
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V- F' 



(5) 



The divergence of the perturbed radiation flux is expressed as 
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where y = 7 + ^Tip^C- in the envelope) and Kp, kj denote the 
logarithmic derivatives of the opacity, while A = (^Ijr') ' the 
(negative) temperature scale height. The unknown temperature 
perturbation T' IT (with its first and second order derivatives 
with respect to r and 0) is then ehminated through the perturbed 
equation of state 



P' 



r 



P 



(7) 



where the x symbols denote as usual the partial logarithmic 
derivatives of the pressure with respect to temperature, den- 
sity and mean molecular weight fi. We have used the fact that 
diffusion may be neglected, so that the Lagrangian perturba- 
tion 6n = 0. Using P' and p', instead of T' and p', gives a 
cleaner numerical solution in the regions with a (not perfectly 
smooth) yu-gradient just outside the convective core. In the stel- 
lar interior the perturbations are almost adiabatic, so that the 
introduced complexities in the energy equation do not mat- 
ter, whereas in the non-adiabatic surface regions the yu-gradient 
vanishes. 
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2.1. Differencing 

The remaining four partial differential equations are then ap- 
proximated by four difference equations on a 2D grid of 3000 
radial and 128 equidistant zones from -QXo njl. The four 
unknown perturbations in the other hemisphere follow from 
symmetry (or anti-symmetry) about the stellar equator, depend- 
ing on whether we force the star with odd or even / - m values. 
The basic 'difference atom' used for the 3-level differencing 
in (r, 9) space is a cross with all four unknowns , ^g, P' /P 
and p'/p defined at the five levels corresponding to the centre 
and the four end points of the cross. This enables the treatment 
of second order derivatives in the energy equation and allows 
for the introduction of turbulent viscosity which, in convec- 
tive regions, is added to the radial- and 9- equation of motion. 
The viscosity is calculated with simple mixing length theory , 
corrected with a reduction term (iGoldreich & Keeleyl (1 1977^ 1 
when the convective timescale is larger than the oscillation 
period. The mixing length is taken as A,„ = min(2 Hp, dmin), 
where c/„„>, = max(|r-r„.|, 0.25Hp) with Hp the pressure scale 
height and r„ the radius of the convective core boundary. The 
stellar model used for the unperturbed (spherically symmetric) 
star was obtained with a version of Eggleton's evolution code 
dPols et al.1119951) . 



2.2. Stability Analysis with forcing 

We use a similar stellar model of a 6Mq star with a core hy- 
drogen abundance of = 0.4 with a radius of 4.3 4 We 
used the same OPAL opacities as IWalker et all (12005!) . With the 
newly developed (see above) 2D oscillation code working with 
a (r, 9) grid it was determined for which complex frequencies 
(T forcing with a specific / and m component of the spherical 
harmonic forcing potential ([T]) yields a maximum (resonant) re- 
sponse. This was accomplished by scanning through a range of 
real forcing frequencies and zooming in on apparent resonant 
frequencies and then optimising the stellar response (signifi- 
cantly) by introducing an imaginary component in the forcing 
frequency. Forcing at the (complex) resonance frequency o-q 
excites (very nearly) the free oscillation mode in which we are 
interested, whereby Im{d-()) is a good estimate of the growth or 
damping rate of the free mode. Forcing was performed with the 
lowest possible I value for a given m value and adopted equa- 
torial (anti)symmetry. We consider only forcing with m=l and 
m=2, so that the forcing has / < 3. The number of nodes ng in 
the range < < ;r of the harmonic forcing functions ^F™ is 
given by ng = / - m. During the zoom-in procedure (close) to a 
resonance the oscillatory response often has a higher spherical 
degree than the forcing, e.g. ng=A for forcing with «e=l-m=0. 
The converged oscillation mode (with complex frequency) may 
then have ng - 2. The excited oscillations for a given m-value 
can in principle have any spherical degree consistent with the 
forcing symmetry (but is limited by the grid resolution). This 
complicates the searching procedure as the search may focus on 
a (stable) higher degree oscillation that occurs in the same fre- 
quency range and is excited through the coupling caused by the 
Coriohs forces. Since it becomes only clear during the zoom-in 
procedure on which kind of mode (spherical degree and stable 
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Fig. 1. The observed (MOST) amplitudes in mmag as a func- 
tion of frequency in the inertial frame / = cr/(2 n) (Walker et 
al. 2005). 



or unstable) the search terminates, this cannot be prevented. 
The Coriolis coupling with higher degrees complicates the sta- 
bility properties of the modes. It appears (see results section 
below) that the resonant overstable modes found are often of 
higher degree than the applied forcing. Some of the modes in 
an instability interval (nmm, nmax) appear stable. This is caused 
by the stabilising coupling with higher degree/order modes. 

2.3. testcase: MOST observations 



Walker et al.l (l2005h give an analysis of the photometric obser- 
vations of the Be star HD 163868 with the Canadian MOST 
satellite in which three clusters of oscillation modes appear, 
see Fig.[T] They modelled stellar oscillations and searched for 
an evolutionary model of the Be star which fits the oscilla- 
tion spectrum: starting with a 6 Mq ZAMS star with X - 0.7 
and Z - 0.02 they arrived at a stellar rotation period of 
^ 0.72 d and evolution away from the ZAMS to a stellar ra- 
dius of R - 4.3/^(7). The opacity was obtained from the OPAL 
1996h . Their analysis (by H.S. and 



tables ftg lesias & RogersI 



U.L, based on a truncated expansion in spherical harmonics) 
indicated that the oscillations with (inertial frame) frequencies 
around 0.02 mHz and 0.04 mHz correspond to prograde, even 
g-modes with \m\ = 1 and \m\ - 2, respectively. They inter- 
pret the oscillations with frequencies near 0.005 mHz as due 
to (retrograde) r-modes with \m\ - 1. Apart from these lat- 
ter modes, they only find pulsational instability for prograde 
modes. The very low frequency oscillations / < 0.005 mHz, 
where f-crliln), could not be identified in their analysis. 

Very recently a new study (iDziembowski et al .l2007h based 
on the traditional approximation (TA) appeared which focusses 
on the oscillation spectrum of the same Be star. Below we dis- 
cuss the results of both earlier analyses and compare them with 
our findings. 
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3. Results 

3.1. Results for symmetric (even) modes 

We adopt a constant angular speed Q = 0.6 Qc (with 
1.70 X lO"'* see'') corre sponding to a rotation period of 0.71 d, 



like Walker et al. I (12005!) and, for some calculations, also Q = 
0.7 Qc . From now on we express all angular frequencies in 
units of Q,.. 
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Fig. 2. The calculated growth rates -Im{&) of overstable even 
oscillation modes versus frequency in the inertial frame for 
Q = 0.6, see Table [T] for mode identification. Plusses: m=l; 
crosses: m=2 



Table 1. Characteristics of the unstable even modes for Q = 0.6 
as shown in Fig. |2] The number of radial nodes (outside the 
convective core) in Im[^r(r, 0)] is given by n^, while rig gives 
the number of nodes in the range < 6* < ;7r of "RelT'ir, 0)]; fi 
and /„ denote the lower and upper frequency of the instability 
interval for the given m- value. 



m 


«s 


fir 






MmHz) 


MmHz) 


1 


2 


21; 25 : 31 


-0.666 


-0.557 


-0.0018 


0.0012 


1 


4 


35 


0.826 






0.0386 


1 


2 


18<" : 40 


0.797 


0.531 


0.0306 


0.0364 


2 


2 


21 : 28; 33 : 34 


-0.675 


-0.543 


0.0135 


0.0178 


2 





14 : 25 


0.418 


0.240 


0.0390 


0.0438 


2 


2 


25 : 31 


0.732 


0.618 


0.0492 


0.0523 



: ng=2 modes with 11^=19,20 are missing. Instead we find two unsta- 
ble ng=4 modes with n,.=40, 42. 



The calculated growth rates versus frequency in the inertial 
frame of overstable even oscillation modes are plotted in Fig.|2] 
for a stellar rotation frequency Q=0.6. Table [T] lists the char- 
acteristics of the unstable even modes found by the 2D code. 
It shows that almost all detected overstable even modes have 
higher angular degree than the applied forcing. The only ex- 
ceptions are the lowest degree (ng = 0) prograde m=2 modes 
with frequencies cr between 0.240-0.418. We find a spectrum 




Fig. 3. Results for g42, an unstable high degree (ng = 4) m=l 
prograde g-mode with frequency f=0.0367 mHz in the ob- 
server's frame. Plotted is the 0-variation of the relative tem- 
perature perturbation T' /T at 10 different radial shells in the 
stellar envelope. The largest amplitude curves correspond to 
the outermost layers. Adopted stellar rotation rate is O = 0.6. 
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Fig. 4. The calculated growth rates - Jm(o-) of overstable even 
oscillation modes versus frequency in the inertial frame for Q - 
0.7, see Table|2]for mode identification. Plusses: m=l; crosses: 
m=2. 



of much higher degree (than the forcing) resonances close to 
the unstable modes listed in tables [T] and |2] Usually these high 
degree/order modes are found to be stable, but not always. In 
Fig.[3]we show the temperature perturbation in the envelope for 
a high degree (ng=4) unstable m=l g-mode, listed in the foot- 
note of table [U The applied forcing has l=m-l, or ng=Q. The 
lowest degree (ng = 0) modes that were found to be overstable 
by Walker et al. (near / ^ 0.02 mHz ) were all stable according 




to our analysis. iDziembowski et al.l (120071) came to the same 
conclusion. Unlike Walker et al., we do not find all overstable 
g-modes to be prograde, although we fin d significantly more 



progra de than retrograde modes unstable. IDziembowski et al 
(l 2007l) . on the other hand, found no significant difference be- 
tween prograde and retrograde modes. However, their analy- 
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Fig. 5. Results for g23, an even m=l retrograde g-mode with ultra low frequency f =0.0003 mHz in the observer's frame. In the 
panel on the left we plot the radial distribution of the local kinetic energy in the mode, while the righthand panel shows the 
variation (with ng - 2) of the relative temperature perturbation T' IT at 10 different radial shells in the stellar envelope. The 
largest amplitude curves correspond to the outermost layers. Adopted stellar rotation rate is Q = 0.7. 



Table 2. Characteristics of the unstable even modes for Q - 0.7 
as shown in Fig.|4] See Table[T]for meaning of symbols. 



m 


ne 








fAmHz) 


MmHz) 


1 


2 


22 : 28; 35 : 36 


-0.703 


-0.562 


-0.0001 


0.0037 


1 


2 


23<" : 40 


0.749 


0.573 


0.0344 


0.0392 


2 


2 


21<" : 34 


-0.727 


-0.578 


0.0182 


0.0222 


2 





15 : 26 


0.387 


0.234 


0.0442 


0.0484 


2 


2 


27(1) . 4^(2) 


0.736 


0.585 


0.0537 


0.0578 



': Some of the modes with = 2 seem to be missing. We do find 
(stable) higher degree modes in these intervals. 
^: the number of radial nodes is not well defined due to lack of resolu- 
tion in the /i-gradient layer 

sis is based on the TA for which substantially more unstable 
modes are found than either by Walker et al. or by the present 
stability analysis. Because the low frequency retrograde modes 
have an extra pair of angular nodes due to the presence of the 
branch of 'buoyant r-modes' the lowest angular order is absent, 
so that there tend to be more prograde than retrograde unstable 
low frequency modes when this is not compensated by higher 
degree instability intervals. The observed cluster of overstable 
modes near f=0.02 mHz are, in our interpretation, high radial 
order retrograde m=2 g-modes. 

We do find unstable low frequency modes with / < 0.005 
mHz that were observed by MOST and could not be explained 
by Walker et al. These modes appear in our analysis as retro- 
grade even m=l g-modes with ng = 2, presumably of higher 
angular degree than sear ched for by Walker et al. Indeed, 
Dziembowski et al.l (l2007b came to a similar interpretation of 
the very low frequency modes. An example of a ultra-low fre- 
quency m=l mode (in observer's frame) with & ^ -Q. can be 
seen in Fig. |5] 

In the left panel of Fig. |5]we plot the local kinetic energy in 
consecutive spherical shells of the grid 



where dr is the radial gridsize. We observe that the kinetic 
energy of the mode peaks in the /i-gradient layer just out- 
side the convective core where the Brunt-Vaisala frequency is 
large and the radial wavelength short. Due to rotational effects 
(Coriolis force) g-modes are confined to an equatorial band 
COS0 < |a"|/(2Q). This effect can be seen in the right panel 
of Fig.|5] where the amplitude in the (evanescent) polar regions 
of the star is almost zero. Although the ^-distribution of the 
temperature perturbation of the retrograde m=l modes (Fig.jS) 
indicates an angular order higher than the forcing, with ng - 2 
(i.e. there are two nodes in the range < < n) and the pos- 
itive and negative parts cancel to some extent. Nevertheless, 
these higher degree modes may still be visible in the photome- 
try even if we observe the star almost equator-on, this has to be 



studied in more detail, see lDziembowski et al.l (120071) . 



LKE(r) = 



f + {y'gf + iy'^f) r^ sinflded^dr 



It appears that several modes (n^ values) inside a given in- 
stability interval are either missing or found stable, while their 
neighbours are again unstable, see tables [1] and |2l This phe- 
nomenon is absent when we apply the TA and use the same 
(now ID) forcing method to analyse the pulsational stability. 
It is caused by the stabilising effect of Coriolis coupling with 
stable higher degree modes, as was also noted by Lee (200l . 
Most of the damping occurs in the /^-gradient layer just outside 
the convective core where the modes have short wavelengths 
and where the local kinetic energy of the modes is maximal. 

Comparing our results with the MOST observations, it can 
be seen in Fig.|2]that the unstable even retrograde m=l modes 
all have inertial frame frequencies / < 0.0012 mHz, while the 
observed distribution in Fig[T]is spread over a larger frequency 
range. This can be improved by adopting a slightly higher ro- 
tation rate. We tried 0=0.7 for which, in the inertial frame, the 
unstable m=l modes are spread out to ^ 0.004 mHz, more like 
what is observed, compare Figs. [T] and H) However, the two 
clusters consisting of unstable retrograde m=2 modes (near f 
^ 0.02 mHz) and prograde m=l modes (near f^ 0.04 mHz) 
then shift to slightly too high frequencies. A rotation rate some- 
where between 0.6 and 0.7 might give the best fit, although the 
observed mode clusters in Fig. [T] remain a bit more extended 
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compared with the calculated spectrum. Our interpretation re- 
quires that the MOST photometry detected only modes with 
even symmetry about the equatorial plane, i.e. HD 163868's in- 
clination angle should be near njl. Cancellation effects would 
then explain that the overstable odd modes (in particular the r- 
modes) remain invisible in a photometric observation. A large 
inclination angle is consistent with the observed large value 
V sin / ^ 250 km/s for this star. 

3.2. Results for anti-symmetric (odd) modes 
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Fig. 7. The calculated growth rates -Jm(o-) of overstable 
modes with odd symmetry versus frequency in the inertial 
frame for Q. = 0.6. See Table [3] for mode identification. 
Crosses: buoyant r-modes with m=l (left) and m=2 (right), 
plusses: g-modes with m=l and asterisks: g-modes with m=2. 



Table 3. Characteristics of the unstable odd modes shown in 
Fig. |7] see Table [T] for meaning of symbols. The retrograde 
modes with the lowest values of \'Re(&)\ are the buoyant r- 
modes plotted in Fig.|7] 



m 


rig 








MrnHz) MmHz) 


1 


1 


18 : 27 


-0.363 


-0.310 


0.0064 


0.0079 


1 


1 


21 : 27 


0.491 


0.424 


0.0295 


0.0277 


2 


3 


35 : 37<" 


-0.709 


-0.690 


0.0133 


0.0138 


2 


1 


14 : 18; 23 : 26 


-0.321 


-0.272 


0.0238 


0.0251 


2 


3 


~ 35 : 37<" 


0.745 


0.716 


0.0526 


0.0518 


2 


1 


17 : 33 


0.679 


0.434 


0.0442 


0.0524 



' : the number of radial nodes is not well defined due to lack of resolu- 
tion in the /j-gradient layer 

Fig Ellists the unstable odd modes found in our calculation, 
see Table [3] for mode identification. For both m=l and m=2 
we found buoyant r-modes (or 'quasi g-modes') destabilised by 
th e /c-mechanism, a res ult which corrobo rates earlier findings 
by ISavonTie (2005) and Townsend (2005) with the traditional 
approximation and bv lLed (120061) . These unstable r-modes are 



the retrograde modes with the smallest \cr\ values in Table[3] In 
Fig. |6]we plot (right panel) the 0- variation of a buoyant r-mode. 
Th e m- 1 (odd) r-mode s near / ~ 0.005 mHz were also found 
by lWalkeretal](l2005b . 

Fig.|6]shows that the high radial order buoyant r-mode ex- 
hibits a 0-variation similar to that of g-modes, although the ra- 
tio of ^r/^e is still significantly smaller than for g-modes. The 
'quasi g-mode' in Fig. |6] shows considerable truncation near 
the pole, almost like a (normal) g mode. For the odd modes we 
found, apart from the buoyant r-modes, only four unstable ret- 
rograde modes (the three asterisks near f ~ 0.013 mHz and the 
plus near f~ 0.005 mHz in Fig.|2ll. 

4. Conclusions 

We interpret the wh ole oscillation spec trum derived from the 
MOST photometry (iWaUcer et al.ll2005h as being modes sym- 
metric with respect to the stellar equator (even modes). It is 
assumed that the modes with odd equatorial symmetry are not 
observed because the Be star is seen almost equator-on, so that 
the light variations of the odd modes cancel to below the detec- 
tion limit. The high value of v sin / ^ 250 km/s suggests indeed 
that the inclination angle / is close to 7r/2. We identify the very 
low frequency (f< 0.005 mHz) oscillation modes in the ob- 
served spectrum (MOST) of HD 163868 as unstable retrograde 
m=l oscillations and the intermediate cluster of modes around 
f=0.02 mHz as retrograd e m=2 oscillations. A similar conclu- 
sion was reached by Dzie mbowski et al.l (l2007h . For Q.-Q.6 the 
highest frequency modes near f=0.04 mHz are identified as 
a combination of prograde m=l and m=2 oscillations. For a 
higher rotation rate of Q=0.7, however, most of the prograde 
m=2 modes s hift ou tside the frequency interval considered by 
Walker et al.l(l2005h . 



The picture that emerges from the analyses by Walker et al 
(W05), Dziembowski et al (D07) and the present study is not 
clear cut in that the detailed results differ in all three studies. 
W05 find a series of unstable prograde even m=l modes that 
are all stable according to D07. We find the lowest angular or- 
der («fl=0) m=l also stable, but find unstable prograde m=l 
modes with ng = 2. W07 find in general much more unsta- 
ble modes than W05 or is found here. This is clearly linked to 
their use of the TA: in our test calculations with the TA (and 
the same input model as for the 2D code) we also find more 
unstable modes, especially retrograde modes, than with the 2D 
code. This seems related to the fact that the TA does not take 
the coupling of different angular degrees by the Coriolis force 
properly into account. Both in W05's analysis and here signif- 
icantly more prograde than retrograde modes are found unsta- 
ble, contrary to the TA findings of D07. 

W05 uses Lee and Saio's method in which the solution 
is approximated by a truncated set of spherical harmonics. 
However, the low frequency modes in rapidly rotating stars can 
deviate significantly from spherical harmonics, so that trunca- 
tion can be problematic. In the 2D code used here no a priori 
assumptions about the 6 behaviour of the solution (apart from 
boundary conditions) are made, i.e. no truncation occurs and 
the solution is limited only by the grid resolution. In this respect 
the present method seems therefore superior to the other two 
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Fig. 6. Result for 123 with f=0.0249 mHz, belonging to a series of unstable m=2 buoyant r-modes. Adopted rotation rate is Q.-Q.6. 



methods and does simulate the observed oscillation spectrum 
not in all details (e.g. the spread of the observed clusters of un- 
stable modes is larger) but the overall result seems quite reason- 
able in view of the simplifications used (spherical, uniformly 
rotating star). The drawback of the implicit 2D method is its 
requirement of fast extended computer memory and patience 
of the user; a single frequency calculation (of which some 20 
are required for a zoom-in on a resonance) takes about 40 min- 
utes on a double Xeon processor node of the compute cluster 
LISA at SARA. 

A point of interest is the visibility of the modes as a func- 
tion of inclination angle, especially those with the higher angu- 



lar degrees, like the progra de m=l modes with no-2. See, e.g. 

et al.l OOOTh for a discussion on the visibility of 



Dziembowski 1 



modes. This would need further exploration, together with a 
study of the effects of differential rotation. 
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